Thermodynamics of Dipolar Chain Systems 
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The thermodynamics of a quantum system of layers containing perpendicularly oriented dipolar 
molecules is studied within an oscillator approximation for both bosonic and fermionic species. The 
system is assumed to be built from chains with one molecule in each layer. We consider the effects 
of the intralayer repulsion and quantum statistical requirements in systems with more than one 
chain. Specifically, we consider the case of two chains and solve the problem analytically within 
the harmonic Hamiltonian approach which is accurate for large dipole moments. The case of three 
chains is calculated numerically. Our findings indicate that thermodynamic observables, such as the 
heat capacity, can be used to probe the signatures of the intralayer interaction between chains. This 
should be relevant for near future experiments on polar molecules with strong dipole moments. 

PACS numbers: 05.65.+b,67.85.-d,05.70.-a 



I. INTRODUCTION 

Rheological electro- and magnetofluids are very flexi- 
ble systems where long-range dipolar interactions lead to 
self-organized filamentation. Since its discovery in fer- 
rofiuids by de Gennes and Pincus this phenomenon 
has been intensely studied in classical systems This 
is typically done using Monte Carlo simulations of hard- 
and soft-sphere dipolar particles Q and comparison can 
subsequently be made to experiments with colloids that 
display many intriguing phases 043 • Important tech- 
nological applications can be found in nanostructuring 
using atomic lithography J8| and in colloidal quantum 
dot photovoltiac devices [0]. However, as the systems 
become increasingly smaller one must start to consider 
also quantum aspects of the structure and the dynamics. 

The explosive experimental development in ultracold 
atoms has lead to a wealth of new possibilities in terms 
of simulating the properties of condensed-matter systems 
with extensive control over both quantum statistics, in- 
teractions, and geometry (Tol - [l^ . Building on this suc- 
cess, a great deal of effort has been put into the cooling 
of Rydberg atoms, permanent magnetic dipole moment 
atoms, homo- and heteronuclear molecules, and samples 
of near-degenerate polar molecules can now be produced 
[l3l - [T5j . These systems provide access to strong dipole- 
dipole forces that are long-range and anisotropic [l6j . 
In particular, in low-dimensional geometries it has been 
shown that the systems are more stable towards two- 
body chemical reaction loss [15[ and to many-body in- 
stabilities [rzi - 

Here we focus on a system consisting of polar molecules 
that are restricted to move in a stack of two-dimensional 
planes. This setup resembles the classical ferrofluids and 
colloidal systems with the major exception that for de- 
generate or near-degenerate ultracold molecules the sys- 
tem is intrinsically quantum and physics beyond the clas- 
sical regime can be addressed. An external field can be 
used to align the dipoles [11], and we assume they are 
always perpendicular to the layer planes. Theoretically, 
a host of interesting few- p^| - |25j and many-body [1^, [23| 



states have been predicted in single- and multilayer sys- 
tems. For very large dipole moments, a triangular crys- 
tal phase has been predicted [1^ that is similar to the 
Wigner crystal of the electron gas [2^ . All these predic- 
tions demonstrate the similarity of the layered dipolar 
case to systems with Coulomb interactions in quantum 
Hall, (S^l graphene, [3l| or semiconductor [12] bilayers and 
to stripe formation in high-temperature superconductor 
materials [33| . 

In the current investigation, we are interested in the so- 
called dipolar chain liquid proposed in Ref. 26 for bosonic 
molecules and later also for fermionic molecules 0, [35| . 
For perpendicularly oriented dipoles, the longest chain 
possible (one molecule in each layer) is always the most 
bound structure [2^. In the case of bosons, it was found 
that the l ong est chain will condense below a critical tem- 
perature [261 ] which depends slightly on whether internal 
vibrational modes of the chain are included |36|. T his has 
been confirmed by a recent Monte Carlo study (37[ , where 
a lattice model was employed. Common to all the stud- 
ies mentioned is that they neglect the repulsive intralayer 
interaction of the polar molecules. Here we include this 
effect in a parametric approach and considered the im- 
plications on thermodynamic quantities such as entropy 
and heat capacity. Thermodynamic variables have been 
addressed experimentally in the case of the unitary Fermi 
gas [111 and for two-dimensional bosonic systems ^39] . In 
particular, density fluctuation measurements can deter- 
mine compressibility [i^j and heat capacity [4l'|. Since 
absorption imaging has been demonstrated also for polar 
molecules [13] , these measurements should be possible for 
such systems as well. 

Below we present results for the thermodynamic be- 
havior of two- and three-chain systems for any number 
of layers. In the case of two chains, we determine the par- 
tition function analytically. Three chains requires, how- 
ever, numerical calculation of the degeneracy of many- 
body states of given total energy. Our formalism is based 
on an effective harmonic approximation for the two-body 
interactions to reduce the Hamiltonian to an exactly solv- 
able form. Naturally, this requires a careful choice of 
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FIG. 1: Schematic view of the setup we consider here for four 
layers with two polar molecules in each layer. All molecules 
have dipole moments that are perpendicular to the layers as 
indicated by the black arrows. The chains are indicated by 
the ellipsoids. The interlayer interaction is denoted Va (at- 
tractive at short distance), while the intralayer interaction, 
Vr, is purely repulsive. In the lower right-hand corner we 
show the choice of coordinate system. 



parameters for the effective harmonic interactions. The 
paper is organized as follows. In Sec. |IT] we discuss the 
Hamiltonian and the various approximations we make. 
As a warm-up, we discuss in Sec.|lll]the case of a single 
chain, before moving onto the case of two- and three- 
chain systems in Sec. IIVI The resulting thermodynamic 
results are presented in Sec. El along with a discussion 
of experimentally relevant systems. In Sec. I VII we sum- 
marize our approach and results, and give an outlook for 
future investigations. 



constant. We will also denote the relative radial distance 
between two molecules by r = y^z^^ + j/^. The interac- 
tion between two molecules in the same layer is found by 
simply setting rt = in Eq. ([T]), i.e. 



Vr{r) = D'-, 



a purely repulsive dipolar term. 

The Hamiltonian for the TV-body system is 



(2) 



H 



2m ^ 



k=l 



dxl 



dyf 



1 



2 ^ 



where {xik,yik) = {xi - Xk,yi ~ yk) 
ordinate of the pair (ifc) with i,k — 



(3) 

is the relative co- 
1, . . . , iV, and the 

factor of i is for double counting, riik is an index that 
denotes the distance (in units of d) between the layers 
containing molecules i and k. For — 0, Vik = K- ^^nd 
for riik > 0, Vik = Va- For simplicity, we assume that 
there is no external confinement in the layer plane itself, 
but this is not essential and can be easily included. Of 
course, it implies that the center of mass motion of the 
total system must be treated separately. 

To solve the A^-body problem, we follow the strat- 
egy outlined in Ref. of approximating the interaction 
terms in the Hamiltonian by quadratic forms for which 
the problem can be analytically solved. Using simpli- 
fied models to describe the behavior of self-assembling 
systems that form complex structures has proven very 
successful 44.] in spite of the diverse origin of the phys- 
ical interactions. For the interlayer potential in Eq. ([T]), 
the harmonic approach has proven very accurate in re- 
producing energy, spatial structure, and wave function 
[23 . |24| . [23, . The real dipole potential is thus replaced 

by 



II. METHOD 

The system we consider is shown schematically in 
Fig. [TJ In general, we have A^ molecules of mass m and 
dipole moment Z? in a layered setup with equidistant lay- 
ers separated by a distance d. The layers are assumed to 
be infinitely thin and allow no interlayer tunneling. This 
will be the case for a deep optical lattice potential. We 
assume that the dipoles are all oriented perpendicular to 
the layers, in which case we have the interlayer potential 



Vaix,y,n) = D^ 



2{ndf 



(a;2 + y^ + {ndffl'^ 



(1) 



where x, y are relative coordinates of two polar molecules 
that are restricted to move in planes located a distance 
nd apart (i.e. adjacent planes have n = 1 and so forth). 
It is convenient to define the dipolar strength, U = , 
and to measure all energies in units of i^j^, and all tem- 



peratures in units of 



id-^kB 



rnd^ ' 

where fc^ is Boltzmann's 



V^''{r,n)=V, 



2{ndY 



- 1 



(4) 



where Vq = Vo{D'^ ,m,d), and the corresponding fre- 
quency (which is denoted a for attractive) is 



2VoiD^,m,d) 
m{nd)-^ 



(5) 



This procedure reduces all interlayer terms to harmonic 
oscillators that are carefully chosen to reproduce the two- 
body energies and the spatial structure of the state (more 
concretely the node of the real potential in Eq. ([IJ is 
reproduced). 

The intralayer repulsion in Eq. ([2]) can also be approx- 
imated by a harmonic oscillator term. However, it must 
be of the opposite sign compared to the interlayer term in 
order to describe a repulsive force, i.e. the term becomes 



2 2 



(6) 
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The problem is now to relate ujr to the properties of 
the repulsive intralayer dipole-dipole potential. This has 
been discussed at length in Ref. [23- 

In the current work, we are interested in the qualitative 
influence that the intralayer terms will have when induc- 
ing an interaction between the chains in the system. We 
therefore use a parametric approach and consider differ- 
ent values of uJr without paying much attention to what 
the physical value may be. However, the systems we con- 
sider are in the self-bound regime, where the attractive 
intcrlayer terms must overcome the intralayer repulsion. 
If we fix uia for given U, then there will be a critical fre- 
quency, LOc, beyond which the A'^-body harmonic Hamil- 
tonian will develop complex solutions. This is the signal 
of instability toward breaking of the chain system into 
smaller systems as discussed in Ref. i25i . 

Here we consider extremal values; uJr = and w,. = 
0.99wc- These choices will help determine what the ef- 
fects of the chain-chain interaction is on thermodynamic 
observables within the harmonic approximation. It is 
then reasonable to expect that the real system will be- 
havior similarly in those limits at least on a qualitative 
level. The major point of the current works is thus to in- 
vestigate how one can determine the role of chain struc- 
tures in the system through thermodynamic measure- 
ments, specifically the changes expected due to quantum 
statistics and due to chain-chain interactions. 



III. SINGLE CHAINS 

We first consider the case of a single chain, which was 
the focus of Refs. [H andlH. Since the geometry is lay- 
ered and no tunneling is allowed, the single chain sys- 
tem is equivalent to a system of N distinguishable par- 
ticles. Solving the harmonic A^-body Hamiltonian yields 
'D{N — 1) (where V is the spatial dimension, V = 2 in the 
case of planes) vibrational degrees of freedom since the 
center of mass separates. We denote the mode frequen- 
cies by uJi, i = 1, . . . ,T>{N — 1). Each of these modes 
will contribute a zero point energy, hiOi. This will be 
part of the overall ground state energy which is not our 
concern here as we focus on the thermal behavior for non- 
zero temperatures and will set the ground state energy 
to zero. 

The partition function in the canonical ensemble, 

0(7V,y,r),is 



(7) 



where gi is the degeneracy of the TV-body state of energy 
El, and /3 = l/(fcBT) is the inverse temperature. For 
the vibrational spectrum of a single chain, each mode 
is non-degenerate and the spectrum is equidistant. The 
partition function is therefore a geometric series and we 
have 



exp(-9,/(2r)) 

i-exp(-e,/r)' 



(8) 



where Qi = tkui/ks- Eq. ^ is the partition function 
of one mode. A state with energy, Ei, involves all the 
modes, and since the energy is in the exponent of Eq. 
we merely need to take the product over all individual 
modes 



exp(-e,/(2T)) 



exp(-e,/T)' 



(9) 



where n is the total number of vibrational degrees of 
freedom in the system, which as mentioned before is 
il = 2{N — 1) when the center of mass has been fac- 
torized. As we will see below, in the case of multiple 
chains, the problem will factorize into a term related to 
single chain behavior, gsop, and a term containing the 
chain-chain modes. 

The center of mass motion can easily be taken into 
account. In the case where there is no external confine- 
ment in the layers, we simply multiply by the well-known 
non-interacting translational partition function 



rnksT 



•D/2 



(10) 



where V is the volume and D is the dimension. The to- 
tal partition function is Q = (Zscpftrans- Ki however, there 
is external confinement in terms of an in-plane quadratic 
field, the number of vibrational degrees of freedom is sim- 
ply VN with T) frequencies that are exactly given by the 
external confinement frequency. The partition function 
can be computed simply from Eq. (j9|). If the system 
isotropic is and the center of mass separates, then Eq. © 
can be simplified to 



/N-l 



q,ep{N,V,T)= n - 



exp(-9,/(2r)) \ 



V 



exp(-e./r) / 



(11) 



Notice that our description of the single chain includes 
all possible internal vibrational modes similar to the ap- 
proach of Ref. m. However, we fix the potential to repro- 
duce physical two-body features of the interlayer dipole- 
dipole interaction, whereas Ref. [sB uses an expansion of 
the dipolar potential around the origin to second order. 
Both approaches yield accurate results for energies of 
longer chains when compared to exact numerics [45j . 



IV. TWO AND THREE CHAINS 

In order to explore the effect of the intralayer repulsion, 
we need to extend to at least two molecules per layer. 
The minimal model would thus be a two-chain system. 
Ideally we would like to consider the thermodynamics 
for any number of chains. However, the requirements of 
quantum statistics makes this an extremely complicated 
problem to solve due to the computational effort involved 
in determining the degeneracy of the many-body states. 
Fortunately, systems of ultracold molecules often have 
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very low density, and in this case we expect the dynamics 
of a few chains to capture the leading correlations in the 
system. The reasons for considering both two and three 
chain systems is that the dipolar interaction is expected 
to induced a tri ang ular crystal in the limit of very strong 
dipole moment [28| . A precursor of this transition should 
be present in a three chain model. For high density, one 
would need a different starting point, typically mean- 
field theory. We concentrate on this low-density limit 
from now on and calculate the properties of two- and 
three-chain systems. 

When dealing with identical particles, one must take 
into account the appropriate symmetry requirements 
of the wave function of those particles. For identical 
fermions, the wave function should be antisymmetric, 
and for identical bosons, symmetric. The ground state 
of the harmonic solution is a symmetric Gaussian wave 
function, which could thus be a good approximation to 
a Bose system. Once we consider excited states for use 
in the partition function, the problem is much tougher, 
especially for many particles, as most excitations will, at 
first glance, appear neither symmetric nor antisymmet- 
ric. However, as we discuss below, in the case of two 
chains, the required symmetrization can be done analyt- 
ically. For three chains one must resort to a numerical 
determination of the partition function [46j . 



A. Two Chains 

For two chains (as shown in Fig. [T|), molecules in dif- 
ferent layers can be treated as distinguishable particles, 
while quantum statistics are important for molecules in 
the same layer. The total partition function can be writ- 
ten as a product of the separable single-chain term of 
Eq. ([9]) and an intralayer quantum statistically restricted 
term 



(12) 



where q^^i, is the boson/fermion intralayer part. 

From diagonalization of the A^-body Hamiltonian, we 
obtain 2W — 1 [i^ frequencies, where W is the number of 
layers, and the center of mass has been separated. Out of 
these frequencies, W — 1 refer to the single chain modes. 
The remaining W frequencies describe the intralayer re- 
pulsion [2^. Each of these W frequencies refer to a 
particular mode in the diagonalizing coordinates. When 
one writes these new coordinates in terms of the original 
laboratory coordinates [43] introduced above, one sees 
that they are in fact antisymmetric in the coordinates 
of the identical molecules in each layer, i.e. they cor- 
respond to traditional relative coordinates, ri^k — i'2,ki 
where fc = 1, . . . , is a layer index. This immediately 
implies that exciting any one of these W intralayer modes 
by an even number of quanta will produce a symmet- 
ric state under exchange of two particles in the layer, 
while an odd number of quanta produces an antisym- 



metric state. In turn, this can be used to construct the 
excitations that are allowed for given quantum statistics. 



1. Bosons 

Bosons require a symmetric wave function under in- 
terchange in each plane. In two dimensions, this can be 
accomplished if all quanta of excitation are either even 
or odd. It can also be done with a combination of even 
and odd numbers, as long as they are the same in each 
dimension for each quantum number, i.e. nix^ niy = odd, 
and n2a;,n2j, = even. The partition function for bosons 
now becomes 

w / 

9^ = n + 1) exp(-2j0fc/T) | , (13) 

where the product runs over all layers, k = 1,...,W^, 
while the sum runs to infinity. This expression can be 
written 



96 



w 

n 



exp(-efc/r) + exp(~39fc/r) 

[1 - exp(-2efe/r)]2 



(14) 



2. Fermions 

For fermions, we require antisymmetry with respect to 
particle exchange. This is achieved by enforcing that the 
total excitation in each direction is odd, i.e. ni^ + niy = 
odd. We have 



w 



If 



which becomes 



9/ 



= 2 



w 



w 

n 



exp(-2efe/T) 



[1 - exp(-20fe/r)]2 



(15) 



(16) 



The prefactor 2^ accounts for all possible interchanges 
between the two directions. For example, for a three 
layer system with two molecules in each layer, there will 
be three intralayer pairs that must be antisymmetrized. 
The lowest fermionic state has an additional degeneracy 
of eight. Using the notation \nixn2xn3x,'n'iyn2yn3y) for 
the quantum state in these three frequencies, the states 
are |111,000), |110,001), |101,010), and |011,100), and 
an additional four from interchange of x and y. Every 
excited state get this additional degeneracy factor of 2^. 



3. Extensions 

The symmetry considerations above can be straight- 
forwardly extended to three dimensions. In addition, an 
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in-plane external confinement potential given by a har- 
monic oscillator does not change anything. The center 
of mass mode will be symmetric in all coordinates so it 
merely modifies the separable part of the partition func- 
tion as discussed above. Thus, it does not modify the 
symmetry properties. 

In the oscillator formalism, a constant can be added to 
the Hamiltonian in order to achieve some desired prop- 
erty of the energy of the system. Here we are interested 
in the behavior of the system as it is thermally excited, 
and not in the ground state energy. We therefore choose 
to set the zero of energy at the ground state level, which 
is of course non-zero in absolute terms due to zero-point 
contributions of all oscillator terms. 



4- Thermodynamic quantities 

Having obtained the partition function, Q, analytically 
in the case of two chains, thermodynamic observables 
can be calculated. In the canonical ensemble, the most 
straightforward is the Helmholtz free energy, F, which is 
the characteristic thermodynamic function of the canoni- 
cal partition function. In the case of distinguishable par- 
ticles, we have 



= ^BT^(^^+ln(l-exp(-e,/r))j, (17) 

where the sum is taken over all vibrational degrees of 
freedom in the system. If the system is isotropic, then 
the free energy is 



w-i 



F = VkBTj2(^+ ln(l - exp(-e,/T)) j . (18) 



The energy, E, is 



E = kuT 



,d\nQ 



dT 

■Dk yVe. , 9,exp(-e,/r) 
^^^^2. (,T+i-exp(-e,/r) 



(19) 



The constant volume heat capacity, Cy, is simply the 
constant volume temperature derivative of the previous 
equation 



c. ^ ill 



dT 

VkB J2 



(20) 



V 

W-l 



e,y exp(-e,/r) 
Tj [1 - exp(-e,/r)p 



The entropy, S, is 



S — kB^^Q + ksT 



d\nQ 
dT 



W-l r 



VkB J2 



ejexp(-ej/T) 

T(l~exp(-e,/T)) 



ln[l - exp(-ej/T)] 



(21) 



The total partition function for indistinguishable par- 
ticles is Q = QscpQ f/b- The partition function can thus be 
separated, and the above formulae can be applied, with 
the appropriate limits on the sums. For bosons, using 
Eq. (fT4|) . the free energy is 



Fb -kBTlnqsep 
w 

-fcsr^ln[cxp(-ej/T) + cxp(-3ej/T)] 



w 



+2fcBT^ln[l -exp(-2ej/T)], 
j 

the entropy is 



(22) 



S 



^sep ~ Ab/T 



W 

E 



kB^j 
T 



exp(-ej/T) -h 3exp(-3ej/T) 



exp(-ej/T) -f exp(-3ej/r) 



4exp(-2ej/T) 



1 - exp(-2ej/r) 
the energy is 

w 



(23) 



E = E. 



Sep 



Y^kBQj 



exp(-ej-/r) + 3exp(-3ej-/r) 
exp(-ej/r) + exp(-3ej/r) 

4exp(-2ej/T) 



1 -exp(-2ej/T) 



and for Cv we have 



Cv, 



^ 4fcBe2exp(-2ej/r) 



Sep 



E 



y2 



(24) 



(25) 



1 



_ [1 + exp(-2ej/T)]2 [1 - exp(-2ej/T)] 
For the fermion case, the free energy is 
Ff = -kBT\nq,,p-WkBT\xi2 



(26) 



w 



^2kBTj2 



T 



ln[l - exp(-2e,/T)] 



The other observables proceed in the same manner as for 
bosons. The factor from the different permutations, 2^, 
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FIG. 2: Helmholtz free energies, F (upper row), and energies, E (lower row), divided by NksT for two bosonic (left column) 
and two fermionic (right column) chains. The different lines correspond to different values of the dipole strength, U, the number 
of layers, W, and the magnitude of the repulsive intralayer frequency, cur, measured in percentage of the critical frequency at 
which the system breaks apart. The nomenclature is {U,W, 



contributes to the entropy, and as a constant factor to 
the free energy, but does not contribute to any observable 
related to any derivative of the partition function. The 
entropy is 



w 



Ss 



2efcexp(-2efc/r) 

T[l - exp(-2efc/r)] 



ln[l-exp(-2efc/T)] +iyfcBln2, (27) 



the energy is 



Ef 



2k, 



w 



29,exp(-2e,/r) ' 
1 - cxp(-2efe/T) 



and finally the heat capacity is 



Cv,sep + 
, W 



rp2 



E 



egexp(-29fc/r) 
[1 - exp(-2efe/T)]^ 



(29) 



B. Three Chains 

The analytical results above for the partition function 
in the two-chain case can unfortunately not be extended 
to a larger number of chains. In this case one needs to 
return to the general partition function Eq. ^ and com- 
pute this numerically. As demonstrated in Ref. !46l, this 
can be done quite effectively in the harmonic approxi- 
mation for both bosonic and fermionic particles. The 
three-chain case has been computed numerically in order 
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FIG. 3: Entropies, S, divided by Nks for systems of two (upper row) and three (lower row) chains of bosonic (left column) 
and fermionic (column) molecules. The nomenclature of the legends is the same as in Fig. [2] 



to obtain g^/^ from Eq. 
obtain Eq. (0). 



m that goes along with qgcp to 



A. Energy 



RESULTS 



We now present the results for both two and three 
chain systems for different values U — 10 and 20, for 
LOr = and w,- — 0.99a;c, and for layer numbers W — 6 
and 12. These values of U are both in the regime where 
the harmonic approximation to Va of Eq. ([T]) is extremely 
accurate. We first discuss Helmholtz free energy and en- 
ergy for two chains, omitting the case of three chains 
which is very similar. Next we present the entropies for 
two and three chains, and finally we show the heat ca- 
pacities. We specialize to also discuss the distinct low 
temperature behavior of the heat capacity where quan- 
tum statistics is most pronounced. 



In Fig. [21 we plot the Helmholtz free energy (up- 
per row) and energy (lower row) for bosons (left col- 
umn) and fermions (right column) for two-chain systems. 
An immediate striking difference between the two up- 
per plots, is the behavior at zero temperature where 
bosons are smooth, whereas fermions have a kink. This 
is a consequence of the ground state degeneracy of the 
fermionic system and appears through the entropy term 
in F = — TS. We will discuss this in more detail when 
we present entropies below. At slightly higher tempera- 
tures, F agrees quite well for bosons and fermions, with 
fermions having lower overall values in all cases. In the 
case of ujr = O.QOwc, the curves are considerably below 
the LUr — case. This is a consequence of the fact that 
the mode frequencies are lowered as 



25|, allow- 



grows 

ing easier thermal excitation of the system. We also note 
that changing the layer number from 6 to 12 induces only 
slight quantitative changes. The energies in the lower 
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row of Fig. [2] show the same tendencies as F, and their 
close similarity attests to the fact that the entropy term 
is causing the differences in F at low temperature. 

B. Entropy 

The entropies of boson (left column) and fermion (right 
column) two (upper row) and three (lower row) chain sys- 
tems are shown in Fig. [31 For two fermionic chains we see 
that the entropy does not go to zero at zero temperature. 
This is caused by the fact that antisymmetry requires at 
least one quantum of excitation. For two-dimensional os- 
cillators the first excited state is degenerate and this im- 
plies finite entropy at zero temperature. From Eq. (|27p. 
the zero temperature contribution is ksW 1112, so the 
plot approaches In 2/2 for many layers (recall that N is 
the total number of particles, while W is the number of 
layers, i.e. W/N = 1/2). For three chains, the ground 
state is non-degenerate in two dimensions (the oscillator 
has a non-degenerate ground state and a twice degener- 
ate first excited state). Thus the entropy goes to zero at 
zero temperature for three chains irrespective of quantum 
statistics. 

We notice that the entropy starts to increase at larger 
temperatures for bosons than for fermions, implying that 
the activation energy gap is larger for bosons. This will 
become more pronounced when we consider heat capac- 
ities below and we postpone the discussion. In general, 
we see that for larger ujr, the entropy in general grows 
for all temperatures as the modes go down in energy as 
discussed above. Again we observe that the quantum 
statistics does not seem to play much role as the left and 
right columns in Fig. [3] are very similar. 

C. Heat Capacity 

The derivatives of the heat capacity makes it very sen- 
sitive to changes in behavior due to various modes that 
become activated in the system. This can be seen for in- 
stance in second order phase transitions where the heat 
capacity will become singular at the critical tempera- 
ture. As we discussed in the introduction, the heat ca- 
pacity is now a measurable quantity in ultracold atomic 
gases [llj. We therefore anticipate that similar measure- 
ments should be possible also for experiments with po- 
lar molecules in the near future. Another good observ- 
able would be the compressibility of the system. The 
detailed balance between attractive and repulsive inter- 
actions should be visible there as well. 

In Fig. m we show a plot similar to Fig. [3] for the heat 
capacities for the same temperature ranges employed in 
Figs. and [3] Overall, we note the high temperature 
behavior of all cases which approaches the equipartition 
value of 2NkB characteristic of two-dimensional oscilla- 
tors. Another general feature is that for larger U, the 
heat capacity decreases due to the larger cost of exci- 



tation in the system. When we compare fermions to 
bosons, we find that fermions in general approach the 
limit slightly slower for both two and three chain sys- 
tems. There is a clear 'shoulder' at low temperature 
which seems to be more pronounced for bosons, but ap- 
pears at lower temperatures for the fermions. Also, going 
from two to three chains seems to make the 'shoulder' 
more pronounced. This feature is clearly associated with 
the larger value of ojr , and also becomes more clear with 
increasing W. This implies that measurements of the 
heat capacity can potentially yield information about the 
role that chain-chain interaction plays in a dipolar chains 
liquid. 

In order to further investigate the interesting features 
at low temperature in the heat capacities, we show a plot 
where the range of temperature is zoomed in around the 
low temperature region in Fig. [SJ Here we see that there 
is striking differences between fermions and bosons, and 
between the two and three chain systems. Not unex- 
pectedly, the effects of intralayer repulsion is most pro- 
nounced at low temperatures. The clear 'shoulder' seen 
for large and three fermions is the one we see at higher 
temperatures for bosons in Fig. 2) This can be under- 
stood from the fact that bosons have larger activation 
energies due to a gap in the low energy spectrum. For 
bosons one cannot generate a totally symmetric state by 
promoting one molecule from the ground to the first ex- 
cited state, since the wave function would then neces- 
sarily become antisymmetric. Thus, the gap for bosons 
is twice the level spacing in the low-energy spectrum, 
whereas for fermions this is allowed and the gap is just 
the level spacing. This is true for the three chain systems 
as seen clearly in the lower row of Fig. [5] For two chains, 
we find that the gap is always the same (twice the level 
spacing) for both kinds of molecules, and the stronger 
increase of Cy for bosons at low temperature is due to a 
larger density of states [i^. Note also the difference for 
three fermions for the two different uJr values, with the 
larger one giving a much smaller activation energy due 
to the lowering of the mode frequencies. 

It is interesting to note that the heat capacity grows 
with W when ujr — 0, while the opposite is true for large 
ojr- This is connected to the fact that in the former case, 
the longer chains have more (degenerate) modes to ther- 
mally excite, whereas in the latter case the repulsion is 
pushing all modes down toward zero but slower for longer 
chains where more attractive terms are present. This 
means that measurements of heat capacities for systems 
of different W can help determine whether the repulsive 
intralayer interaction is important or not for given system 
density. 



D. Relation to Current Experiments 

For comparison to experiment it is convenient to re- 
call our units for dipolar strength and for temperature. 
Here we use that 1 Debye equals 3.336 • 10~'^°Cm in SI 
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FIG. 4: Heat capacity, Cv, divided by Nks for systems of two (upper row) and three (lower row) chains of bosonic (left 
column) and fermionic (column) molecules. The nomenclature of the legends is the same as in Fig. (2] 



units. Also we must remember that the dipole-dipole 
interaction in SI units has a factor of where en is 

the permittivity of vacuum. The dimensionless dipolar 
strength can be written as 

where u is the atomic mass unit. Similarly, the unit of 
temperature that we use can be expressed in the form 

The experiments most relevant to our setu p a re the 
ones being conducted at the moment at JILA [1^. The 
rotational and vibrational ground state molecules used 
are '*'^K^^Rb which has a maximum dipole moment of 
D = 0.566 Debye. The experiment reported in Ref. [HI 
has studied chemical properties of a sample confined to a 
quasi-two-dimensional stack of about 23 layers and found 



profound effects due to the optical lattice potential on the 
loss rates. The dipoles are perpendicular to the planes 
and the applied field has a magnitude that induces a 
dipole moment of D = 0.158 Debye. With a lattice 
spacing of 532 nm, this amounts to U = 1.15 for maxi- 
mum dipole moment, while for the experiment cited it is 
roughly U = 0.09. Our results are based on the strong- 
coupling harmonic approximation, and U must therefore 
be large. The JILA experiments are therefore not in the 
regime of validity of our approach. 

To get a stronger dipole moment, NaK molecules con- 
sisting of a Sodium and a Potassium atom can be used. 
They are predicted to have a maximum moment of about 
D ^ 2.7 Debye [l^. This has the advantage that several 
isotopes of Potassium are available for ultracold atomic 
experiments of both bosonic and fermionic character. Ex- 
periments on this system are currently under way [48| . 
Assuming that d = 532 nm, one could in principle reach 
values of L/ ^ 13 which is within the region where we ex- 
pect our approximations to hold. Furthermore, the range 
of temperature shown in Figs. above corresponds to 
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FIG. 5: Same as Fig. [4] but for a smaller range of temperatures around zero. 



27 nK in this case. We thus expect that intralayer effects 
could be studied for system temperatures below about 
270 nK. Reaching this scale should not be difficult ex- 
perimentally. 

VI. SUMMARY AND OUTLOOK 

The physics of dipolar chain systems in multilayered 
geometries holds great potential for realization of exotic 
few- and many-body physics that is hard to access in 
other setups. Here we have considered the case of chains 
of dipolar molecules with dipole moment oriented perpen- 
dicular to the layer planes and studied the thermodynam- 
ics for both bosonic and fermionic molecules. In order to 
render the problem tractable, we approximate the real 
Hamiltonian by a harmonic form with suitably chosen 
parameters that reproduce realistic two-body properties 
of dipolar systems. 

In contrast to previous studies of the dipolar chain liq- 
uid, we include the effect of interactions between the 
chains by introducing a repulsive harmonic term that 



parametrizes this effect. We have also explicitly taken 
care of the quantum statistical effects of having more 
than one molecule per layer. In the simplest non-trivial 
case, there are two chains in the system and for that we 
can determine the partition function and in turn all ther- 
modynamic properties analytically within the harmonic 
approach. In the case of three chains, the degeneracies 
of the states could no longer be determined analytically 
and we resort to numerics. 

Our findings indicate that heat capacity measurements 
for the system can probe the relative influence of excita- 
tion within single chains and the inter-chain dynamics. 
While we have only considered the two and three chain 
cases here, for a low density system we expect these to 
be the leading correlations in the system. A viable way 
to study our predictions could be through measurements 
at different densities. Starting from very low density, we 
would expect that the single chain modes are more impor- 
tant and will dominate the thermodynamics. However, 
as the density is increased, one should start to see effects 
of the inter-chain modes. In the case where the repulsive 
intralayer interaction is very strong we find a softening of 
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the modes of the system and thus an increase in the heat 
capacity that could be observable. Another observable 
that we expect will be strongly influenced by inter-chain 
dynamics is the compressibility. In the limit of very large 
density, we expect our approach to become less accurate 
and mean- field theory is probably a better starting point. 

The melting of classical crystals including those with 
dipolar interactions have been extensively studied (49| . 
and recently the corresponding quantum system in the 
case of a single layer with perpendicular dipoles was ad- 
dressed [1^. The interest is focused on the ability of 
dipolar systems to form a crystal akin to the Wigner 
crystal of the two-dimensional electron gas [l^. The re- 
sults presented here suggest that in the multilayer case 
there can be an interesting competition between inter- 
layer attractive interactions and the in-plane repulsion 
that could modify the crystallization. This is an inter- 
esting topic for future studies. 



Another question that arises for a system consisting of 
dipolar chains is whether one can develop an effective de- 
scription that considers the chains themselves as the basic 
constituents. This is conceivable in the limit of low tem- 
peratures where the occupancy of higher modes is small. 
The quantum statistics of the chain would then depend 
on whether it has an even (bosonic) or odd (fermionic) 
number of molecules. At first glance, the former would 
constitute a (repulsive) Bose gas [5^, while the latter 
would correspond to a Fermi liquid-like state [IH. The 
effective interaction between these entities could help dis- 
cern potential instabilities in the system, such as the ro- 
ton instability for bosons or the density-wave instability 
for fermions. In the limit where the number of layers 
is large the two types of instability should merge. It 
would be interesting to try to come up with an analyt- 
ical expression for the dynamics in this limit. Here the 
harmonic approximation could be very useful. 
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